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. Abstract 

Network coding and cooperative communication have received considerable attention from the 

research community recently in order to mitigate the adverse effects of fading in wireless transmissions 

. and at the same time to achieve high throughput and better spectral efficiency. In this work, we design 

and analyze deterministic and random network coding schemes for a cooperative communication setup 

with multiple sources and destinations. We show that our schemes outperform conventional cooperation 
^ ! 

\ in terms of the diversity-multiplexing tradeoff (DMT). Specifically, it achieves the full-diversity order at 
00 . 

. the expense of a slightly reduced multiplexing rate. We establish the link between the parity-check matrix 

• ■ for a {N + M, M, + 1) systematic MDS code and the network coding coefficients in a cooperative 

' communication system of N source-destination pairs and M relays. We present two ways to generate 



o 



the network coding matrix: using the Cauchy matrices and the Vandermonde matrices, and estabUsh 
that they both offer the maximum diversity order. 



X 

' Index Terms 

Cooperative communication, network coding, outage probability, diversity-multiplexing tradeoff. 



I. Introduction 

Channel fading is one significant cause of performance degradation in wireless networks. 
In order to combat fading, diversity techniques that operate in time, frequency or space are 

H. Topakkaya and Z. Wang are with the Department of Electrical and Computer Engineering, Iowa State University, Ames, 
lA, 50011, USA (e-mail: hakan, zhengdao@iastate.edu). 



April 20, 2010 



DRAFT 



2 

commonly employed. The basic idea is to send the signals that carry same information through 
different paths, allowing the receiver to obtain multiple independently faded replicas of the data 
symbols. Cooperative diversity tries to exploit spatial diversity using a collection of distributed 
antennas belonging to different terminals, hence creating a virtual array rather than using physical 
arrays. 

In [[U Ahlswede et al. introduced network coding to achieve the max-flow rate for single- 
source multicast that could be impossible to achieve by simply routing the data. Since then, 
network coding has been recognized as a useful technique in increasing the throughput of a 
wired/wireless network. The basic idea of network coding is that an intermediate node does not 
simply route the information but instead combines several input packets from its neighbors with 
its own packets and then forwards it to the next hop. However, since network coding is devised 
at the network layer, error-free communication from the physical and medium-access layer is 
usually assumed, which is a simplifying assumption for wireless communications. 

Efforts have also been made to apply network coding to the physical layer, e.g. in [fTOll . [|23l . 
[[24|. Towards that goal, cooperative schemes have been proposed that make use of network 
coding in a cooperative communication setup, and studies have been conducted to determine 
whether network coding provides any advantages over existing cooperative communication 
techniques [H, 0, [El, [23, |l20l-|l22l. 

In [[T2|. a network-coded cooperation (NCC) was proposed and its performance was quantified 
using the diversity-multiplexing tradeoff analysis which was originally proposed for multiple 
antenna systems in f25\. NCC was shown to outperform conventional cooperation (CC) schemes 
which includes space-time coded protocols [81 and selection relaying 1*311: It requires less 
bandwidth, and yield similar or reduced system outage probability while achieving the same 
diversity order. However, these results are based on an optimistic assumption that any destination 
node should receive the packets that are not intended for it without any error so that the intended 
packet can be recovered from the xor'ed packet sent by the relay. When this assumption is 
removed the scheme can no longer achieve the full diversity order of M + 1, where M is the 
number of cooperating relays, but only a reduced diversity order of 2. 

In this paper, we propose a network coded cooperation schemes for source-destination pairs 
assisted with M relays. The proposed scheme allows the relays to apply network coding on the 
data it has received from its neighbors using the coefficients from the parity-check matrix of a 
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MDS code. A closed form expression for the outage probability is derived. We also obtain the 
diversity-multiplexing tradeoff performance of the proposed scheme under two different traffic 
network models: multicast and unicast. Specifically, it achieves a maximum diversity order M+1 
at the expense of a slightly reduced multiplexing rate. We also propose two different network 
coding approaches: deterministic and random. We show that our scheme outperforms NCC and 
CC in terms of probability of outage. 

The rest of the paper is organized as follows. Section|ll]discusses the system model, description 
of the proposed scheme. In Section IIII-Al performance analysis is established using DMT and 
the main result is presented. Section |IV] discusses the network code design. In Section |Vl we 
discuss unicast, random network coding and selection relaying. In Section |Vl] the performance 
of the proposed scheme is compared in terms of DMT and average outage probability with the 
existing schemes in the literature. Section IVIII contains the conclusions. 

II. System Model 

A. General System Description 

The network studied in the paper is composed of source-destination pairs denoted as 
(si, (ii), . . . , (sat, cLm), and M relays denoted as ri, . . . , in a single-cell where all the nodes 
can hear the transmissions of each other as shown in Fig. [IJ We assume that each packet is 
composed of L bits: hi = [&i,i,6j,2, • • • ,&i,L]- We divide bi into smaller blocks of equal length 
/ and represent the k^'^ block [bi^ki+i, &j,fci+2, • • • , k e {1, . . . , K} a finite-field element 

Oi,k & where g = 2' and K = L/l. Therefore, each packet is represented as a -ft"-tuple 
©i = [di,i,0i,2, ■ ■ ■ ,Oi,K] G F^^^'; see e.g., [HI, H. Dividing each packet into small blocks 
enables us to work with a smaller field size which in return significantly reduces the complexity 
of the arithmetic operations. This is to be contrasted to the scheme in [|3 where the field size is 
taken to be g = 2^. We will give a lower bound on the field size in Sec. |lVl We consider two 
different transmission scenarios. In the first scenario, each source node Sj is trying to transmit 
the data packet 6j to all the destinations di,i = 1, . . . , N which is known as the multicast 
scenario. In the second scenario, each source node Sj is trying to transmit the data packet 6i 
to only destination and we will refer to this scenario as the unicast scenario. All the nodes 
are assumed to be equipped with half-duplex (i.e. cannot transmit and receive at the same time) 
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single-antennas. Each data packet 6j is error control coded and modulated, and transmitted in 
T time slots. 

The channel between any pair of nodes is assumed to be frequency flat fading with additive 
white Gaussian noise (AWGN). Let x,; G C denote the transmitted symbols from node i and 
yj E C the received symbols at node j. The additive noise Zi ~ CN(0, 1) has independent and 
identically distributed (i.i.d.) circularly symmetric entries. Let hij E C denote the instantaneous 
channel realization. We assume that the channel coefficient hij remains constant during the 
transmission time of a packet. Then, the channel within one block can be written as 

yj{t) = ,/phijx,{t) + Ziit), t = 1, 2, . . . , T. (1) 

where p is the average received SNR at the destination. All the transmissions are made with 
equal power. In the above equation, the transmitter could be any of the sources or relays, the 
receiver could be any of the relays or destinations, as long as the transmitter and receiver are 
different (i.e., not the same relay). The channel coefficient hij between any two nodes is modeled 
as i.i.d. with zero-mean, circularly symmetric complex Gaussian random variables with common 
variance Therefore, is exponentially distributed with parameter /3 Vz, j . 

A total of NL bits are transmitted by all sources in (A^ + M)T channel uses, therefore the 
system rate is i? = A^L/[(A^ + M)T] bits per channel use (BPCU). The transmission rate Rq for 
one source or one relay per one packet is fixed, identical, and equal to Rq = L/T = R(N+M) /N 
BPCU. 

The instantaneous mutual information of the channel model in ([U) when i.i.d Gaussian input 
is used is given by: 

/(X,;y,)=log(l + |/i„|V)- (2) 

where Xj and Yj denote the transmitted symbol by node i and received symbol by node j. 
We assume that powerful enough channel codes can be applied within each packet such that if 
I{Xi; Yj) > Rq, the packet can be decoded correctly. In case errors occur, we assume they can 
be detected. This can be realized through cyclic redundancy check (CRC) code or other parity 
check codes. When /(X^; Yj) < Rq, we say that the channel hij is in outage. Otherwise, we say 
that the channel hij is operational. Define r = [2[^(^+*^^]/^ — l]/p. Since is exponentially 
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distributed, the outage probability for the channel in ([T]) is given by: 

Po = Pr{I{Xf, Y,) < Ro) = Pr{\h,J^ < r) = 1 - exp(-/3r) - /3r, (3) 
where a(r) = 6(r) if lim^^o[a(r)/fe(r)] = 1. 

B. Network Coded Cooperation 

Our transmission scheme consists of two stages; see Fig. [21 In the first stage, direct 
transmissions from the sources to the destinations take place in orthogonal time slots. Thanks 
to the broadcast nature of the wireless medium, all the destinations and the relays overhear the 
transmissions. At the end of the first stage, each relay tries to decode all N packets. Here one 
of the two strategies is possible: 

1) Strategy £/: If a relay can successfully decode all the packets, then it participates in 
the second stage. Otherwise, it remains silent. In the second stage, the participating 
relays perform network coding. Specifically, relay i will transmit the linear combination 

2) Strategy I^: If a relay can successfully decode at least one packet, then it participates in 
the second stage. Specifically, if relay i was able to decode the packets correctly from the 
sources in the set Si where Si C {1,...,A^}, then it will transmit the linear combination 

Unless otherwise specified, we study the first case when the Strategy £/ is used until Sec. IVl 
Strategy ^ will be discussed in Sec. |Vl 

C. Deterministic and Random Network Coding 

We will consider two network coding schemes for the user cooperation: random coding and 
deterministic coding. In the random coding approach, which we will refer to as Random Network 
Coded Cooperation (RNCC), relay draws aij randomly from the finite field F^. After the 
random coefficients are drawn, a new packet is created by making a linear combination of the 
source data packets using the a^/s. In the deterministic approach which will be referred to as 
Deterministic Network Coded Cooperation (DNCC), the coefficients a^/s are predetermined and 
they are designed in a way to maximize the probability that the received linear combinations 
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are actually decodable at the destination. We will discuss the problem of how to choose these 
predetermined coefficients in detail in Sec. HVl 

In order to express the overall transmitted signal, we define the following matrix: 



A 



1 ... 



«M,1 



T 



(4) 



... 1 ai^N ■ ■ ■ OiM,N 

where (■)^ denotes transpose. Also define the N x K finite field vector corresponding to the 
original source packets as 6 = [9f , Gf^, . . . , Q^]^. Using matrices A and 6, we can express the 
potential transmitted signals by all the sources and M relays, in that order, as 11 = AQ where 
n G Note that 11 represents the potential transmitted signals, since due to severe 

fading some of the channels might be in outage and therefore only a subset of packets can be 
successfully decoded by some relays. Under Strategy £/, such relays will not participate in the 
second stage and the rows of A corresponding to these relays can be considered to be deleted. 
Under Strategy however, only the coefficients in A that correspond to the unsuccessful 
packets would be zero, as opposed to a whole row being deleted. Note that, from the destination 
di's perspective, some of the channels might also be in outage. We denote the corresponding 
submatrix of A for destination di by Ai which satisfies Hi = AiQ where Ilj denotes all correctly 
decoded packets at destination di. 



III. Performance Analysis 

A. Diversity-Multiplexing Tradeoff 

As mentioned in the introduction, we will investigate the performance of the proposed scheme 
via diversity-multiplexing tradeoff (DMT). DMT is accepted as a useful performance analysis 
tool in cooperative systems [[3]|, [[T2l|. For completeness, we give the formal definitions as in 
[|25l . Let -Pe(p) denote packet error probability of user i at SNR p. Define = minj P^, 
i = 1, . . . , N, then a scheme is said to achieve spatial multiplexing gain r and diversity gain 
d if the data rate is limp_^oo R{p) / log(p) = r, and the minimum error probability satisfies 
limp^oolog(Pe(p))/log(p) = -d. 
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B. Main Result 

Next we define a new parameter which plays a key role in the derivation of the outage 
probability and hence the achieved diversity order. For any integer i e [1, min(m, n)], we define 
the F-rank, Fi(C), of a m x n matrix C as an integer 7 such that 1) any collection of 7 rows 
of C is at least rank i, and 2) there exists a collection of 7 — 1 rows of C that has rank i — 1. 
Next, we derive the DMT of the system as a function of Tn{A). 

Theorem 1. The diversity-multiplexing tradeoff of DNCC with N source-destination pairs and 
M intermediate relay nodes which choose their linear combination coefficients from the matrix 
A for multicast using Strategy is given by: 



d{r)^{N + M-{rr,{A)-l)) 



N + M ' 



Proof: 

1 ) Multicast: In the multicast problem, the necessary and sufficient condition for destination 
di to recover ©j is rank{Ai) — N. To analyze the outage probability, we define the following 
events: Ei = {rank(Ai) < N}, and E^^ = {A has at most F;v(^) - 1 rows}. Notice that, 
Ei C -B"^ by the first condition in the definition of F-rank. By the second condition in the 
definition of F-rank, there exist a collection of rows of A that are rank N — 1. Let Ai denote 
a {Tn{A) — 1) X N submatrix of A that consists of such rows. Let denote the event that 
m relays fail to receive all the Gi's correctly. Define El""" = {Fq n {A, = A,}}. It follows that 
0OW Notice that the probability that any relay can successfully decode all packets in 

the first stage is P{S) = UtiP^^sA^-^Y) > ^0) = n,=iexp(-^T) = exp(-7V^T). As a 
result, 

P(F^) = (^) P(5)^-'"(l - P{S)r. (6) 

Having N direct transmission from the sources and M — m transmissions from the relays, 
each destination can potentially receive and decode N + M — m packets. Let E{N + M — m,l) 
denote the event that I out of N + M — m channels were operational: 



(7) 
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where Pq is given by ©. Since El""" C Ei C E^^, using ^ and (|7]) we have: 

M rjv(A)-l 

P{E,)<P{Ef^) = J2PiFm)- P{E{N + M-m,l)) (8) 

m=0 1=0 

and 

PiE,) > PiEl""") = P(Fo)Po^^^"^^'^^^^"'^(l - Po)^'^^'^^^' (9) 

In ([8]), the first summation stands for the probability of the event that m of the relays fail to 
receive all Bj's correctly, leaving us with only M — m relays which will participate in the second 
stage. In total N + M — m transmissions will be made. The destination di may not be able to 
recover all 0i's, if only Vn{A) — 1 or less number of transmissions are successful. 
Notice that, as p — )■ oo, r — )■ 0. We need to find the following limit: 

iTo r^+M-(r.(A)-i) - (10) 

We consider the individual terms in the summations one-by-one and find the term with the 
smallest order of r. Observe that lim^_^o (1 - P{.S)) = N(3 and P{Fm) = Kmr"^ where Km is: 



r->-o \ m 

fk\ pk-l 



Similarly, P{E{k,l)) = K^^it^ ' where K^^i = The smallest order r term happens 

when / is equal to Tn{A) — 1. Hence, we have: 

P(PJ P{E{N + M-m,l)) _ , 

I™ ^N+M-m-irr,{A)-l) " ^«^-^N+M-'m,^^{A)-l - ^up, (12) 

and 

Similarly, we can show that P(Pf"') = K,o^r^+*^-(r^(^)^i) where Ki^y, = /3^+A^-(riv(A)-i)_ 
Now, choosing the fixed rate to be P = r log p and substituting into (fT2l) . we obtain: 

P(Pi) = irp(T^-i)(^+*^-(riv(A)-i)) (^4) 

where -ft'iou; K < Kup, which is the desired result. 
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2) Unicast: In unicast we have a different problem: given the received packets Fj at destination 
di, we would like to recover only 0j from the set of linear equations = AjO. The error can 
only happen when the direct link is in outage. Notice that this implies that Ai does not contain 
Cj (the i'th row of the N x N identity matrix Inxn)- In this case, a necessary and sufficient 
condition for 6j to be recoverable is that G span{Ai), where span(Ai) is the row-space of A^. 
Here, we make another rank definition that will be useful for the proof of the unicast scenario. 
We define the Aj-rank, Aj(C) of a m x n matrix C as an integer A such that 1) any collection 
of A rows of C spans a space that contains Cj but 2) there exists a collection of A — 1 rows of 
C which does not span a space that contains Cj. Next, we derive the DMT of the system as a 
function of Ai(A). 



Lemma 1. DMT of DNCC for unicast for the i^^ destination is 

d, {r) = {N + M-{K{A)-l)) 



N + M 

1 r 

N 



r e ( 0, — 1 . (15) 

^ ' N + M ' ^ ^ 



Proof: Here we define the following relevant events. Ei = {cj span{Ai)}, Ej^^ = {A^ 
has at most Ai{A) — 1 rows}. Notice that, E^ C -Ej"^ by the first condition in the definition of 
Aj-rank. By the second condition in the definition of Aj-rank, there exist a collection of Aj(A) — 1 
rows of Ai that does not span e^. Let A^ denote a (Aj(y4) — 1) x submatrix of A^ that consists 
of such rows. Keeping the definition of define e/"^ = {Fq fl {Ai = Ai}}. It follows that 
e/""" C Ei. Since E]""" C Ei C Ei'"'', using © and © we have: 

AI Ai(A)-l 

P(E,) < P(E,"n = Po 5^ i'li^J ■ Yl PiE{N-l + M-m,l)) (16) 

m=0 /=0 

and 

P{E,) > PiEt"") = P(Po)Po'^+*'"^^'^^^"'^(l - Po)^'(^)-' (17) 



where the first Pq in (1161) accounts for the outage of the direct link between Sj and dj. The limits 
in the second summation in (fT6l) is due to the fact that the destination di may not be able to 
recover all Bj's, if only Aj(A) — 1 or less number of transmissions are successful. The rest of 
the proof can be completed by showing that the diversity orders of both the upper and the lower 
bound are equal to (fT5l) as in the proof of the Multicast scenario. ■ 
Notice that from the definition of F-rank, we have r7v(^) = maXjAj(A). Since the error 
probability Pe is defined to be the minimum of individual error probabilities, we have d{r) = 
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mmidi{r). After substituting di{r) in (fTSl ). we obtain the desired result in ([5]). ■ 
Corollary 1. The maximum diversity is achieved if and only ifVN^A) = N. 

Proof: The result follows immediately from dS]). ■ 

IV. Design of the Linear Network Coding Matrix 

In this section, we try to design a network coding matrix A that can yield the maximum 
diversity order. Notice that, by definition we have Tpf(A) > N. Therefore, it is clear from 
Corollary [T] that we need to pick an A that satisfies Tm{A) = N. Before going into the discussion 
on the design of the matrix A, we would like to give another important rank definition that will 
be used in the design of A. 

The row Kruskal-rank [fTSl . ifTSl of A, denoted by k(A), is the number r such that every set 
of r rows of A is linearly independent, but there exist one set of r + 1 rows that are linearly 
dependent. 

Lemma 2. k{A) = N -^TnIA) = N. 

Proof: We prove T^iA) = N ^ f^iA) = N; the other case is straightforward. When 
Tn{A) = N, from the definition of F-rank any collection of rows of A is at least rank A^. But 
since rank{A) < N, we have the first condition for the Kruskal-rank. Also since rank{A) < N, 
any + 1 rows will be linearly dependent. ■ 
The minimum Hamming distance dmin between any two codewords for a (n, k) error-correcting 
code is upper bounded by the Singleton bound as dmin < n — k + 1. The codes that achieve this 
bound are called maximum distance separable (MDS) codes [fTTTl . The following result relates 
the column Kruskal-rank of the parity-check matrix if of a linear block code to its minimum 
distance dmin' dmin = f^{H) + 1. This follows from the following theorem in IfTTl by realizing 
that dmin = n — k + 1 for an [n, k, d] MDS code ^. 

Theorem 2. ( fTl\ p. 318]) 'if is MDS if and only if (iff) every n — k columns of H are linearly 
independent. 

The transpose of the parity check matrix of a systematic (A^ + M, M, A^ + 1) MDS code 
can be used as an encoding matrix A for our DNCC scheme to minimize the total number of 
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packets necessary at the destinations for decoding the source packets. If such an A is used, 
then each destination needs and only needs N packets (from the sources and relays) for correct 
decoding. Note that depending on the sizes N and M, finding a {N + M, M, N + 1) MDS code 
may or may not be possible in a given finite field Fg IfTTI . |fT6l . 



A. Network Code Designs from RS Codes 

Reed-Solomon (RS) Codes are MDS codes. There are two ways of constructing an RS Code: 
either using the Vandermonde matrices [HU or using the Cauchy matrices fl9\ Sec. 4.3]. Because 
of the special structure of A in (HI), we will be working on systematic RS codes. 

1) Construction based on Cauchy Matrices: The systematic generator matrix for the RS(n, k) 
code has the form G = [I\C] where / is the identity matrix of order k and C is a k x (n — k) 
matrix [fT9ll and G satisfies k{G^) = n — k. C is known as Cauchy matrix and is given by: 



UiVj 



0>z>A;-l, 0>j>n-k-l. 



Xi + Vj 

where Ui,Vj,Xi and yj are elements of GF{2"^) and are defined as: 



Xi 



/3 



n—l—i 



0> i> k-1. 



%• = /3"-'-'-^ 0>j>n-A;-l, 
1 



0> i> k-1. 



(18) 

(19) 
(20) 

(21) 
(22) 



0>i>fc-l 

where (3 is the primitive element for ¥g. Therefore, choosing n = N + M and k = M, we 
construct the network code A = [I\a] by choosing aij = Cij which gives fi(A) = N. 

2) Construction based on Vandermonde Matrices: The Vandermonde matrices are defined 
from a vector of m distinct generating elements {ti, . . . ,tm} of Fg as: 



1 u ti 



t 



n—1 



It t^ 



j.n— 1 



(23) 
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The determinant for the square Vandermonde matrix is given by det{Vnxn) = Y[i<i<j<n('^j 

and Vn is nonsingular iff all the t/s are distinct. To construct A for given N and M, we do the 

following: 

1) Choose a suitable ¥g with q = 2^ > N + M. 

2) Choose N + M distinct elements ti, ^2, • • • ,'tN+M of Fg. 

3) Construct the Vandermonde matrix Vatxa^ from ti, . . . , t^v and Vmxn from tAr+i, . . . , t^+M- 

4) Then = iyMxNV^lj^)i^j and A = 

Note that the generating elements that are needed in the construction of RS codes from 
Vandermonde matrices requires an extra property that they should be the consecutive powers of 
the primitive element (3 G Fg, i.e. . . . , for a t-error correcting RS code. We do not 

need or impose this requirement. 

Lemma 3. Let T ^W^^"- be an invertible matrix. Then, for any H G F™^" is the same as that 
ofHT: Ti{H) = Vi{HT) and k{H) = k{HT), V 1 < i < min(m,n). 

Proof: Pick i G (1, min(m, n)) arbitrary rows from H and denote the resulting n x n matrix 
by H' . We need to show that if the rows of H' are linearly dependent or linearly independent 
then so are the rows of H'T. But notice that, since T is full-rank we have xH' = -v^ xH'T = 
where a; G F^ and F^ is the iV-tuples in Fg. ■ 

Picking H = [Vn^nIVmxn] and T = Vj^^j^ and using Lemma[3]we have rjv(G') = k{G) = N. 
Note that since we need N + M distinct elements of the finite field F^, it is enough to have 
q > N + M. Next, we give an example for the case when N = 2 and M = 2. 

Example: Consider a [n, k, d] = [4, 2, 3] MDS code <^ over F4 = {0, 1, a, /3 = = a + 1} 
with symbol representations as {0 = (0,0), 1 = (0, l),a = (1,0), /3 = (1, 1)}. Constructing the 
Vandermonde matrices from the set {0, and multiplying with the inverse, we have 



Let K = 1 and bi = [bi^i, 64,2], ©i = [^j.i], ^ = {1)2}. If the above encoder matrix is used, relays 
will transmit the linear combinations /36'i 1 + 062^1,061^1 + (362^1 respectively. For example for 
the first relay, this operation will be performed using regular addition in F2 as (61 2 + &2,i + 
^2,2, &i,i + &i,2 + &2,i) and for the second relay (61^1 + 61^2 + ^2,2, &i,i + &2,i + ^2,2) which can be 




(24) 
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put in the matrix form as: 








1 


1 


1 


\ 


/ 




\ 




1 


1 


1 









bi,2 






1 


1 





1 






b2,i 




V 


1 





1 


1 


/ 


V 


b2,2 


J 



(25) 



Notice that since k{A) = 2 and the destination will be able to recover 6i, 62 when at least 
two of the transmissions are successful. It is important to emphasize that unlike MDS code 
construction that we gave earlier, here we do not have any restrictions on the code size for any 
given and M, we can find a large enough finite field ¥q satisfying q = 2\l > N + M. 

V. Discussions and Further Improvements 

In the previous sections, we have established the DMT of DNCC for a given network coding 
matrix A. Later, we designed the network coding matrix A to have the property that Tn{A) = 
k{A) = N. Here in this section, we look at the performance of DNCC under the Strategy 
^ and the case when only some of the relays (which are selected according to their channel 
qualities) are allowed to transmit. We also investigate the performance of the case when the 
linear combination coefficients are chosen randomly. 

A. Strategy SS 

Decode-and-forward schemes suffer from performance loss when the source-relay channel is 
in outage. If a multi-source scenario is considered the performance loss becomes even more 
severe. Therefore, the assumption that the relay has to decode all the packets in order to be able 
to cooperate may be too restrictive for such schemes. We could relax this assumption and assume 
that the relays will participate cooperation even though they have not been able to decode all 
the packets. 

Denote the outage event under Strategy SS by Ef . Notice that under Strategy not only 
the M — m relays as in ^ but also the rest of the m relays contributes to the decoding at the 
destinations. Clearly the probability of not being able to solve the linear system of equations will 
decrease and hence the performance will get better, i.e. P{Ef) < P{Ei). Taking Tn(A) = N, 
the probability of El""" becomes P(Ef"') = P{Fo)P^^+\l - Po)^-\ We have PiE^""") < 
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P(Ef) < P{Ei). But using a similar analysis as in the proof of Theorem [H it can be shown 
that even though Strategy ^ offers lower packet error rate, the DMT is the same as that of 
Strategy . That is, even though Strategy improves the packet error rate performance, the 
DMT remains unchanged. 

B. RNCC: 

In RNCC the linear combination coefficients «j,/s are chosen randomly from a finite field 
¥q. Similar to the deterministic case, destination rfj cannot recover when the submatrix Ai is 
rank deficient, i.e. Ei = rank{Ai) < N. However, unlike the deterministic case there are two 
possible reasons to have a rank deficient in the random case: one is due to fading and the 
other is due to the choice of the random coefficients ai/s. The former happens when at most 

— 1 channels are operational resulting in an Ai matrix that has at most — 1 rows. Notice 
that no matter what a^^'s are chosen A^ will be rank deficient and hence the linear system 
of equations cannot be solved. Therefore, we define this event to be deterministic error event: 
]j>det _ Yi^^ at most — 1 rows} in the proof of Theorem [T] and probability of this event is 
given by P(Ef *) = ^^^^^ P(F„) ■ ^^^^ P{E{N + M - m, I)). Notice that * C Ei. On the 
other hand, due to the random choice, relays may choose linearly dependent coefficients which 
will result in an Ai matrix such that rank{Ai) < N. Denote this event by E'-"-"-. But by the 
Corollary [H this event will result in the outage events that have diversity order less than M+ 1. 
Therefore, the key idea of this proof is to isolate such events, and show that the probability of 
such events can be bounded by the field size. 

P{Ei) = P{Ei\Ef''*)P{Ef^') + P{Ei\El'''')P{El'''') 

= P{Ef') + P{Ei\El'''')P{El'''') < P{Ef') + P(Ei|E™") (26) 

Next, we present the lemma that upper bounds the term P(£'j|£'[""). 

Lemma 4. The probability that any N x N square submatrix A'- of Ai is rank deficient is upper 
bounded by, 

PiEilEl""") < — (27) 
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Proof: The proof consists of an application of the Schwartz-Zippel Lemma of a carefully 
chosen error event as in the proof of Theorem 2 in [|6||. We skip the proof due to space limitations. 

■ 

Now, using Lemma |4] we see that if the field size is large enough, the error event will be 
dominated by the event Ei with high probability: 

lim P{Ei) < P{Ef^') + lim — = P(Ef *). (28) 

q—^oo (J— >oo q 

Note that, although the limit is taken asymptotically with g — )■ oo, it is enough to have 
q = p^^+^ The rest of the proof is the same with the above proof for DNCC. Below we 
summarize all the above proved results with the following theorem: 

Theorem 3. DNCC with M intermediate relay nodes which choose their linear combination 
coefficients from the rows of A that satisfies V]^{A) = N and N source nodes achieves the DMT 
in both the multicast and unicast scenario and under both strategies ^ and SS: 



d{r) = (M + 1) 



\ (N + M) ■ 



RNCC achieves the same DMT as in (1291 ) with probability at least 1 — where q is the field 
size. 



C. Selection Relaying: 

We can also consider the case where not all of the M relays transmit, but only K selected 
relays transmit. The same selection rule based on the instantaneous wireless channel conditions 
can be adapted as in [[T2l|. Define 

hi = iam{\hsiirif , l^ridiJ^, • • • , l^sjvnl^, l^ridivT} (30) 

where hij is the channel coefficient between node i and node j. Then, select the K relays that 
maximizes hi, namely first choose r with the rule: r = argmax^, /ij and continue the same 
process of choosing the maximum in the beginning of each relay transmission. Note that this 
selection mechanism can be implemented using a distributed protocol at the network layer as in 
O: relays choose a timer that is inversely proportional to the quality of their channels. Relays 
assess the quality of their channels from the RTS-CTS packets that were transmitted by the 
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source and destination nodes respectiveljQ. No CSI is required at the physical layer. Next, we 
give the diversity-multiplexing performance of this scheme. 

Theorem 4. DNCC scheme with the selection of the best K relay nodes out of M and N source 
nodes, rjv(^) = in the multicast scenario achieves the DMT: 

<iW = (A- + l)(l-ii^r), ^e(o,-^) (31) 

if K < N — 1, and otherwise achieves the DMT: 

( N + K \ f N \ 

d{r) = {N + M{K-{N-l)))[l-^^ry r e (32) 

Proof: Let r = argmax/ij where hi is as in (l30l) . The cdf for \hjr\^ (or \hrj\'^) where j can 

be a source (or a destination) node was derived in lfT2l as: 

Mr M ^ 

m=l I jjtm J 

where (3m,k^ are the parameters of the exponential random variables associated with the 
corresponding channels between node m and node k, and 13m = Yl!k=i [l^rn,k + /3fc,m] , ^ = 
{1, . . . , M}. Taking f3m,k = (3 and using exponential expansion a high-SNR approximation for 
(l33l) can be shown to be equal to: F(r) = (2A^/3)*^"^/3r^''. Now, the probability that any relay 
can successfully decode all N packets in the first stage is P{S) = Y[k=i ^ (\^kr\^ > r) = 
Y[k=i ~ -^(''")) = (1 ~ ^i'T'))^ ■ Similar to the definition of let denote the event that 
k out of K relays fail to receive all the packets: P{Fk) = (^)P(5)^^'=(1 - P{S))K Using 
similar techniques as in the proof of Theorem [H it can be shown that Pr{Ek) = Kit^'^'' where 

Also let Es^t{k) denote the event that s channels out of N source channels and t 
channels out oi K — k relay channels were operational. Then we have P(Es^t{k)) = 
{^JP^-'{1 - Poy{^;'')F{T)^-''-\l - F(r))*. Notice that, since rjv(A) = we have 
E, = EfP. Therefore, we have P{E^) = Ef=o^(^fe) " E{s;tUt=o} P(Es,t{k)). It can be 



^Forward and backward channels between the relays and the destinations are assumed to be the same due to reciprocity 
theorem (13 \. 
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shown that P{Es^t{k)) = K2T^^-'^^^^^-^~^^^^ , where = max(x, 0) and K2 = 

^A^^)^x-fc^)^iv-s+(i^-/c-t)Af(2iv)(M-i)(x-fc-t) hence P{E,) = K^K2T'ikHN-s)++(K^k-t)+M _ 

We need to find out mmi^st\s+t={o....,N-i}}{{N — s)^ + {K — t)^M). We need to consider 
two different cases: K < N — 1 and K > N — 1. For the first case, choosing t = K 
ands = N - 1 - t = N - 1 - K achieves the minimum: (iV - s)+ + {K - tyu = 
{N - {N - 1 - K) + {K - K)M = K + 1. And for the second case choosing t = N - I 
and s = achieves the minimum: (iV - s)+ + {K - tyu = N + {K - {N - 1))M. Now, rest 
of the proof can be completed using similar techniques as in the proof of Theorem [TJ 

■ 

In [l3l|, for a single source single destination setup it was proved that instead of transmitting 
from all the M relays, if a selection is performed and only the relay with the best channel 
coefficient transmits then the BER at the destination enjoys a M-fold diversity gain. Inspired by 
this idea, the authors in [12] proposed a Network-coded cooperation (NCC) scheme for N s-d 
pairs and M relays where only the "best" relay selected according to (l30l) . XOR'es all the source 
packets and transmits to the destination (Fig. [2] (b)). However, the M-fold diversity order can 
only be achieved when an unrealistic assumption is made. The assumption is that the destination 
has to be able to decode all the other source packets successfully. When no such assumption 
is made, no gain from the selection process is obtained and only a diversity order of one is 
achieved. The significance of our result in Theorem |4] is that if enough number of relays could 
be used (K > N — 1), we can achieve N + M (K — (N — 1)) diversity order. 

VI. Comparison vv^ith Other Schemes 

A. DMT Comparison 

In this section, we would like to compare diversity-multiplexing tradeoff of the previously 
proposed schemes in the literature. The closest scheme in the literature is the NCC scheme 
considered in [[T2|. In NCC instead of all the relays, only one relay transmits which results in 
total of + 1 time slots. Using fewer time-slots NCC achieves a better spectral efficiency than 
the proposed scheme here. However, NCC can only provide a fixed diversity order of two, while 
the the proposed scheme achieves the full-diversity order of M + 1. 

In the following, for comparison we include the DMT performance of NCC and that of conven- 
tional cooperation (CC) which includes space-time coded protocols (HI and selection relaying [3J. 
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The diversity-multiplexing tradeoff of NCC is given by d (r) = 2 (1 - r{N + l)/N) , r G 
('-'' A^)- "^^^ DMT of the decode and forward strategy with M intermediate relay nodes is 
given by [8]: d (r) = (M + 1) (1 - 2r) , re (0, 0.5) . 

To show the advantage of the proposed schemes, we present diversity-multiplexing tradeoff 
of the existing schemes and the proposed schemes in Fig. [3l As can be seen from the figure, 
both of the proposed schemes and CC provide a full diversity order of M + 1 when r — )• 0. 
However, the proposed schemes can provide a higher diversity gain than CC when the spectral 
efficiency increases. 

B. System Outage Probability 

Here, we compare the system outage probability of the proposed schemes with the other 
schemes. The system outage occurs when any di is unable to decode reliably: = 1 — 
YliLi (1 ~ -P(-^i))- We compare Pg with the system outage probabilities (30), (43) derived in 
lfT2l . In Fig. |4]for a network consisting of two source-destination pairs with various number of 
relays is compared. As can be seen from the figure, the proposed method clearly outperforms 
NCC by achieving the full diversity of M + 1 as compared to NCC's fixed diversity order of 
2. The same performance is observed in Fig. [5] where the number of source-destination pairs is 
increased to three. 

C. Monte-Carlo Simulation 

Here, we compare these schemes with the existing schemes via Monte-Carlo simulations for 
various number of relays. In the simulations, only channel conditions are considered to isolate 
the diversity benefits of the scheme. We generate an (N + M) x and an A^ x A^ matrix 
that contains the channel coefficients for each destination and each relay, respectively. Then, we 
decide that the transmission is successful for any link if the instantaneous channel condition 
is large enough to be able to support the given data rate and we update the same size linear 
coefficient matrices accordingly. After all the transmissions take place, we perform Gaussian 
elimination on the updated linear coefficient matrices to conclude whether each destination di 
was able to recover the source packet 6j or not. The channel coefficient variances are chosen 
to be equal to one and Rq = 1 BPCU. Please note that we considered the average outage error 
probability which is found by dividing the total number of errors occurred by the number of 
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source nodes instead of the system error probability. Please note that, in all the figures only 
the unicast scenario is adapted since CC cannot be implemented in a multicast scenario. We 
compare the proposed schemes in Fig. [6] and Fig. |7] with NCC and CC. As can be seen from 
the figures, the proposed schemes are able to provide the M + 1 diversity order and outperform 
the other schemes. 

Theorem [3] claims that the performance loss incurred due to the assumption under Strategy 
is not in terms of diversity gain but it is in terms of coding gain. This is validated through 
simulations as shown in Fig. [8] and Fig. |9] using Monte-Carlo simulations which are set up in 
the same way as in Fig. [6] and Fig. Ul 

VII. Conclusions 

In this paper, we have proposed two network coded cooperation schemes for N source- 
destination pairs assisted with M relays. We studied two different traffic network models: 
multicast and unicast. The proposed schemes allow the relays to apply network coding on the 
data it has received from its neighbors. We allow the relays to linearly combine the packets with 
coefficients either drawn randomly from a finite field or drawn from a linearly independent set 
of coefficients. We established the diversity-multiplexing tradeoff performance of the proposed 
schemes, and showed its advantage over the existing schemes. Specifically, it achieves the full- 
diversity order M + 1 at the expense of a slightly reduced multiplexing rate. The novelty 
in the proposed scheme is that we establish the link between the parity-check matrix for a 
(A^ + M, M, N + 1) systematic MDS code and the coefficients to perform network coding in a 
cooperative communication scenario consisting of source-destination pairs and M relays. We 
presented two ways to generate the network coding matrix: using the Cauchy matrices and the 
Vandermonde matrices. 
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Fig. 1. System model: N source-destination pairs and M relays 
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Fig. 2. time-division allocation for the different schemes compared: (a) CC (b) NCC (c) DNCC, RNCC 
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Fig. 3. DMT comparison of various schemes 
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Fig. 4. System outage probability comparison. (N=2 M=l, ...3) 
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Fig. 5. System outage probability comparison. (N=3, M=l, ...3) 




Fig. 6. Average outage probability per destination, (N=2, M=l, ...3) 
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Fig. 7. Average outage probability per destination, (N=3, M=l, ...3) 




Fig. 8. Average outage probability per destination, (N=2 M=l, ... 3) 
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Fig. 9. Average outage probability per destination, (N=3 M=l, ...3) 



